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Abstract 

RULISON  was  an  underground  nuclear  explosion  detonated  to  stim¬ 
ulate  gas  production  as  part  of  the  PLOWSHARE  program.  The  RULISON 
device  was  detonated  at  21:00:00.12  on  10  September  1969.  It  was 
detonated  in  Western  Colorado,  which  represented  a  new  source 
region,  and  was  emplaced  at  a  depth  nearly  double  that  of  previous 
American  explosions. 

Our  analysis  includes  study  of  seismograms  from  18  LRSM,  2 
VELA  observatories,  LASA,  and  8  WWSS  stations.  When  travel- times 
are  compared  to  the  Herrin  68  tables,  interval  velocities  along  an 
easterly  profile  generally  agree  with  previously  published  findings. 
The  most  probable  magnitude  estimates  based  on  amplitude  data  are: 

JUjj  (Gutenberg)  =  4.62  and  (Gutenberg)  =  3.99. 

Proposed  identification  criteria,  developed  mostly  from  NTS 
explosions  and  Western  United  States  earthquakes  were  applied  to 
RULISON  and  were  compared  when  possible  to  a  previous  chemical  ex¬ 
plosion  at  CLIMAX,  Colorado,  and  to  a  selected  Colorado  earthquake. 
The  criteria  include  location,  depth  of  focus,  complexity,  M  vs 

nij^,  energy  relationships  among  phases,  first  motion,  and  radLtion 
patterns . 

Computed  locations  displaced  the  event  from  the  true  epicenter 
from  1  to  24  km  depending  upon  the  number  of  recording  stations 
used  and  the  source  depth  restrictions  placed  on  the  computations. 
Cepstral  analyses  produced  a  depth  estimate  of  about  2.5  km  for 
RULISON  with  reasonable  velocity  assumptions  although  the  depth 
phase  pP  was  not. visible.  For  the  9  August  1967  earthquake  cep¬ 
stral  analysis  produced  a  depth  estimate  of  5.0km.  At  present, 
however,  we  lack  sufficient  experience  with  earthquakes  to  say  that 
cepstral  peaks  of  the  magnitude  observed  might  occur  for  many  of 
them.  Thus,  strictly  speaking,  the  depth  determination  cannot  be 
used  as  a  discriminant.  Complexities  are  not  applicable  to  RULISON. 


Mg  versus  and  energy  relationships  among  phases  place  RULISON  in 
the  explosion  population;  and  show  that  there  is  less  shear  energy 
present  for  RULISON  than  for  the  9  August  1967  Colorado  earthquake. 
No  rarefactional  first  motion  was  observed  on  RULISON  or  earthquake 
recordings.  The  long-period  LQ/LR  ratios  can  be  explained  by 
surface-wave  radiation  patterns  from  a  compressional  source  accom¬ 
panied  by  tectonic  strain  release  of  relative  strength  0.6. 
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Introduction 

The  underground  nuclear  detonation  RULISON  was  a  part  of  the 
PLOWSHARE  program  to  investigate  gas  stimulation.  However,  by 
virtue  of  its  location,  shot  medium,  and  depth  it  provided  an  im¬ 
portant  set  of  new  data  relevant  to  underground  nuclear  test  detec¬ 
tion  and  identification.  Its  setting  was  in  Western  Colorado 
(39°21'21"N,  107‘’56'53"W),  a  completely  new  source  region  for  in¬ 
vestigation  and  comparison  with  the  other  source  regions  where 
underground  nuclear  detonations  have  occurred.  The  depth  (8,431 
feet)  was  approximately  a  factor  of  two  greater  than  any  previous 
underground  nuclear  explosion  in  the  United  States,  and  the  shot 
medium  (Mesa  verde  Sandstone)  is  different  from  previous  shot  media 
with  the  exception  of  GASBUGGY.  Thus,  the  RULISON  event  provided 
an  excellent  opportunity  to  test  the  effectiveness  of  the  proposed 
diagnostic  criteria  which  were  developed  largely  on  the  basis  of 
NTS  experience. 

In  this  study  we  have  made  an  evaluation  of  a  selected  set  of 
RULISON  seismic  data  and  examined  the  applicability  of  each  of  the 
current  purposed  diagnostic  criteria  for  identifying  underground 
nuclear  explosions.  A  more  complete  data  analyses  by  the  USC§GS  is 
in  preparation.  In  addition  there  was  a  large  chemical  explosion 

(CLIMAX)  in  the  same  general  region  available  for  comparison  with 
RULISON. 

The  basic  observational  data  are  presented  in  a  single  section 
and  saved  permanently  on  magnetic  tape  at  the  Seismic  Data  Labora¬ 
tory  in  order  to  provide  a  readily  accessible  data  base  for  each 
aspect  of  this  investigation  as  well  as  for  future  work.  These 
data  are  then  used  in  the  detailed  evaluation  of  location,  source 
depth,  and  other  diagnostic  parameters  for  the  RULISON  event  as 
compared  to  other  underground  nuclear  detonations  and  earthquakes. 
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BASIC  OBSERVATIONS  AND  MEASUREMENTS 


Tabulated  dAta  '  \ 

-  ,  ■  1  ] 

Epicenter  data  for  the  nuclear  event  RULISON  are  given  in 
Table  I.  Data  for  28  stations  are  presented  in  Table  II.  These 
include  travel-times  for  observed  short-period  phases:  P,  Pg,  and 
Lg.  Maximum  zero  to  peak  amplitudes  (A/T)  and  corresponding 
periods  for  short-period  and  long-period  phases;  and  body-wave  (m. 
and  m^)  magnitudes  and  surface-wave  (M^  andlM^)  magnitudes  at  those 
stations  for  which  they  could  be  computed.  ' 

^  The  stations  listed  ihclude  18  LRSM  stations,  2  VELA  observa- 
torifes  and  8  WWSS  stations.  The  stations  are  listed  in  order  of 
increasing  distance  from'  the  RULISON  epicenter.  The  geographic 
coordinates  and  elevations  for  these  stations  are  given  in  Appendix 

1.  Figure  1,  shows  the  event  location  and  the  LRSM  and  VELA  sta¬ 
tion  network. 


Instrumentation 

The  LRSM,  VELA  and  WWSS  stations  consist  of  three  component 
short-period  and  long-period  instnljmentation.  Pertinent  recording 
inforijiation  for  these  stations  is  givenUn  Appendix  2.  Relative 
magnification  curves  for  all  these  networks  are  given  in '.Appendix  3. 

‘  \ 

Reduced  travel- tim^s  \  ' 

!  1 

Reduced  travel-times  are  shown  in  Figure  2  relative  to  Herrin 
68  times.  We  also  show  the  Pn  and  P  velocities  which  fit  the  data. 
iThese  indicated  velocities  are  used  in  conjunction  with  relative 
Pn  and  P  amplitudes  for  magnitude  estimates  according  to  Evernden 
(1967)  and  are  discussed  in  the  magnitude  section.  , 


Amplitudes 


Figures  3  through  7  are  plots  of  maximuM  amplitude  (zero  to 
peak  A/T)  of  Pn,  P,  Pg,  Lg,  LQ,  and  LR.  For  the  shjort-perioid  Pn 
ind  P  phases,  the  maximum  peak-to-peak  excursion  within  the  first 
few  cycles  of  motion  of  the  phase  was . measured,  halved,  and 


divided  by  the  period  of  the  measured  amplitude  cycle  to  obtain 
zero-to-peak  A/T.  For  both  the  short  and  long-period  Pg,  Lg,  LQ, 
and  LR  phases  the  maximum  peak-to-peak  excursion  visible  on  th^ 
record  was  measured,  halved,  and  divided  by  the  period  ito  again 
obtain  a  zero-to-peak  A/T  value.  Further,  Lg  and  LQ  measurements 
were  made  from  the  records  of  the',  horizontal  instrument  most  nearly 
normal  to  the  direction  of  the  propagating  wave. 

The  expected  attenuation  r^tes  for  Pn  and  P  are  indicated  with 
thd  amplitudes  on  Figure  3.  These  are  discussed  in  more  detail  in 
the  Magnitude  section.  Attenuation  rates  of  R"^  fitted  visually 
for  Pg  and  Lg  data  are  shown  in  Figures  4  and  5.  For  LQ  we  visually 
fitted  a  slope  of  -1.66  through  the  data  points.  Figure  6.  In 
Figure  7,  LR  amplitudes,  we  show  two  slopes  of  -1.16  and  -1.66  for 
distances  less  than  15“  and  a  slope  of  -1.66  for  distances  greater 
tl'jan  15*.  These  slopes  are  based  on  the  surface  wave  magnitude 
estimates  and  are  discussed  in  more  detail  in  the  Magnitude  section. 

In  general,  the  Pg,  Lg,  and  LQ  amplitude  data  conform  to  the 
expected  attjsnuation  rates. 

\ 

Magnitudes  ' 

Body-wave 'magnitudes  are  estimated  in  this  report  in  two  ways 
for  stations  in  the  distance  range  2*-100°. 

mjj  =  logA/T  + 

where 

A  =  maximum  zero  to  peak  ground  motion  (millimicrons)  of  the 
Pn  or  P  phase  within  the  first  t|iree  or  four  cycles, 

'  T  =  apparent  period  of  mai^imum  ground  motion  (seconds),  ^ 

B  =  distance  correction  factor. 

I 

The  distance  correction  factors  from  16*  to  100*  were  taken, from 
Gutenberg  and  Richter  (1956).  The  VELA  Seismological  Center  (VSC) 


1 


B  (Gutenberg  and  Richter  1956) 


projected  this  curve  back  to  10®  and  assumed  an  inverse-cube  rela¬ 
tionship  for  the  amplitude  decay  of  the  Pn  wave  from  10®  to  2®. 
These  factors  are  listed  in  Appendix  4. 

Gutenberg  and  Richter’s  definition  of  body  wave  magnitude 
at  teleseismic  distances  is  accepted  by  practically  all  seismologi- 
cal  organizations  as  the  standard;  however,  Evernden  (1967)  shows 
that  for  distance  less  than  20®  it  is  necessary  to  adjust  the 
distance  correction  factor  in  their  formulation  for  the  Nevada  Test 
Site.  Evernden' s  formulas  used  in  this  report  are  as  follows: 

m?  9  =  -7.55  +  1.21  (log  (A/T)  +  3.04  logA) 

nig  g  =  -3.27  +  log  (A/T)  +  2  logA 


where 


A  =  distance  in  kilometers. 

Application  of  these  formulas  are  based  on  the  appropriate  veloc- 
ites  and  relative  amplitudes  of  P  (Figures  2  and  3).  In  addition 
in  Figure  3  we  show  what  the  amplitudes  should  be,  for  the  various 
magnitude  determinations  relative  to  mj^  (Gutenberg,  A>16®)  »  4.62. 

It  should  be  noted  that  the  Evernden  formulation  of  m,,  was 

not  used  in  these  estimates.  Even  though  the  travel- times  for  LAO, 

TUC,  and  BP-CL  show  an  apparent  8.1  km/sec  velocity,  these  stations 

do  not  correspond  to  the  appropriate  region  or  path  indicated  by 

Evernden  for  this  correction  and;  therefore,  m^  g  was  used  for  the 

magnitude  estimates.  CR2NB  could  be  either  an  m^  ,  or  m,,  ^  accord- 

^•ao«5 

ing  to  Figure  2;  however,  the  amplitude  shown  in  Figure  3  indicates 
that  an  mg  g  would  provide  a  better  estimate  with  respect  to  the 
teleseismic  magnitude  estimate.  The  signals  for  LON  and  COR  appear 
to  travel  at  the  7.9  km/sec  velocity  but  the  amplitude  at  COR  is 
large  and  consequently  the  mg  g  is  applied  at  this  site.  SJ-TX  and 
WQ-IL  show  an  apparent  velocity  of  about  10.5  km/sec;  however,  the 
amplitudes  indicate  that  mg  g  provides  the  better  magnitude  esti¬ 
mates.  We  show  all  magnitude  estimates  in  Table  I  and  Figure  8. 
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On  the  basis  of  Figures  3  and  8,  the  best  body-wave  magnitude  esti 
mates  are  (Gutenberg  A>16“)  m^^  =  4.62  and  the  adjusted  mj^  =  4.59. 
However,  due  to  the  uncertainties,  notei  above,  for  applying 
Evernden’s  corrections,  the  most  acceptable  magnitude  estimates  is 
mjj  (Gutenberg  A>16“)  =  4.62. 

Calculation  of  the  surface-wave  magnitude  for  RULISON  is  by 
the  method  of  Gutenberg  (1945)  for  all  distances.  In  addition  we 
computed  for  distances  less  than  15°  by  a  method  of  Von  Seggern 
(1969). 

Mg  “  1*818  +  C  +  D,  Gutenberg  (1945), 


where 


Apj  =  0.5  peak  to  peak  amplitude  in  microns  at  T  =  20  seconds 
for  the  horizontal  radial  component  of  Rayleigh  wave. 

1.656  logA  =  distance  correction  factor  with  A  measured  in 
degrees.  This  correction  is  limited  to  A  between  15°  and  130°. 

C  *  Site  correction  factor. 

D  =  Depth  of  source,  azimuthal  correction,  etc. 

For  this  study  we  set  C  and  D  equal  to  0  and  use  the  following  re¬ 
lation  adopted  by  Geotech  (1964); 

~  log  *  1.66  logA  -  0.18  where  (A^)  =  peak  to  peak 

amplitude  in  my  and  T  =  corresponding  period  in  seconds  for  the 
vertical  component  of  the  Rayleigh  wave  and  A  =  distance  in  degrees 

These  two  formulas  are  identical  at  T  =  20  seconds.  The 
Geotech  formula  does  not  consider  ellipticity  (A^^/A^);  however,  for 
periods  of  15  to  17  seconds  and  an  ellipticity  of  0.8  the  variance 
of  Mg  (Geotech)  to  (Gutenberg)  is  +  0.03  magnitude  units.  Since 
small  magnitude  events  traversing  continental  paths  have  their 
maximum  measurable  amplitudes  spanning  this  range  of  periods,  the 
formula  is  compatible  with  Gutenberg's. 
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It  Is  important  to  note  that  this  magnitude  estimate  gives  a 

magnitude  of  0.18  less  than  the  USC8GS  and  Bashams  estimates.  Iheir 
formulas  are  as  follows: 

USC5GS,  Mg  =  log  (A/T)  +  1.66  logA  +  3.3 


where 


A  -  distance  in  degrees 


and 


A/T  =  amplitude,  zero  to  peak  in  p/sec; 
Basham,  =  log (A/T)  +  1.66  logA  +0.3 
where 


A  =  distance  in  degrees 


and 


A/T  -  amplitude,  zero  to  peak  in  mp/sec. 


Since  Gutenberg's  formulation  is  for  distances  greater  than 

IS  ,  we  use  a  modified  distance  correction  factor  for  surface  waves 

for  those  less  than  1S«  (Von  Seggern,  1969)  derived  by  the  use  of 

ayleigh  wave  amplitude  measurements  from  29  Nevada  Test  Site  ex¬ 
plosions. 

The  formulation  is  as  follows: 

Mg  =  log(A/T)  +  1.16  logA  +  0.74. 


These  estimates  are  tabulated  in  Table  I  and  shown  in  Figure  9. 
Further,  in  Figure  7  we  show  the  Rayleigh  amplitudes  and  the  ex¬ 
pected  amplitudes  for  Mg^  =  4.18,  Mg  =  3.99  (A>15»)  and  Mg  -  3.84 
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for  all  stations;  however,  the  scatter  of  the  data  makes  it  diffi¬ 
cult  to  determine  which  slope  best  fits  the  data  for  distances  less 
than  15®.  From  Figure  9,  it  is  clear  that  formulation  over 
corrects  and  the  Gutenberg  formulation  applied  to  distances  less 
than  15®  under  estimates  the  surface  wave  magnitude  relative  to 
those  estimated  at  teleseismic  distances.  Therefore,  we  conclude 
that  the  best  magnitude  estimate  is  (Gutenberg  A>15®)  =  3.99. 

The  following  table  summarizes  the  magnitude  estimates  for 
RULISON.  The  word  "Adjusted"  preceding  the  magnitude  symbol  indi¬ 
cates  the  nearer  station  magnitudes  are  corrected  using  either 
Evernden  or  Von  Seggern's  formulas  and  are  included  in  the  average. 
The  mjj  or  for  "all  distances"  includes  in  the  average,  the  mag¬ 
nitudes  determined  by  the  use  of  VSC  distance  correction  factors 
for  mj^  and  Gutenbergs  formula  for  for  distances  less  than  15®. 


(All  Distances) 

(Gutenberg,  A>16®) 

Adjusted 

(All  Distances) 

(Gutenberg,  A>15®) 

Adjusted 


mj^  =  4.90  ^  0.62 

mj^  =  4.62  ^  0.36 

mj^  =  4.59  +  0.35 

Mg  =  3.84  +  0.39 

Mg  =  3.99  +  0.44 

M  =  4.09  +  0.36 
s  — 


for  27  stations, 
for  11  stations, 
for  27  stations, 
for  26  stations, 
for  12  stations, 
for  26  stations. 


RULISON  seismograms 

Various  profiles  consisting  of  LRSM  and  VELA  station  seismo¬ 
grams  for  RULISON  are  shown  in  Appendix  8,  Most  of  the  indicated 
profiles  are  not  profiles  in  terms  of  a  continental  linear  array 
type  due  to  the  paucity  of  stations  and  lack  of  azimuthal  alignment. 
However,  the  East  profile  consists  of  six  LRSM  stations  (CR2NB, 

BY- 10,  WQ-IL,  GZ-OH,  AS-PA,  and  PJ-PA)  which  are  about  equally 
spaced  between  957  and  2757  kilometers  and  on  an  80®  azimuth  from 
RULISON.  Included  also,  is  HN-ME  even  though  this  station  is  not 
exactly  on  the  profile  (Figure  1). 

In  Figure  10  we  show  the  seismograms  for  the  East  profile 
which  indicate  a  Pn  velocity  of  7.9  km/sec  and  later  P  velocities 
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o£  8.4,  10.4,  and  13.1  km/sec.  These  velocities  appear  reasonable 
when  compared  to  the  NE  profile  from  BILBY,  Archambeau  et  al  (1969) 
and  data  from  the  Lake  Superior  study  by  Glover  and  Alexander 
(1970).  Further  the  8.4  km/sec  lies  between  the  8.32  and  8.52  km/sec 
for  the  ESE  and  NNE  profiles  from  GASBUGGY  (Rasmussen  and  Lande, 
1968).  Thus  the  velocities  seen  along  this  profile  are  as  expected. 

Summary 

Fcr  the  RULISON  event  we  have  obtained  the  travel- time  and 
amplitude  data  for  18  LRSM  stations,  2  VELA  observatories,  LASA, 
and  8  WWSS  stations.  These  are  saved  on  magnetic  tape  at  the 
Seismic  Data  Laboratory.  In  addition,  the  seismograms  for  the  LRSM 
and  VELA  observatory  stations  were  digitized  and  saved  for  subse¬ 
quent  analysis. 

We  conclude  that  the  magnitude  estimates  most  appropriate  1( 
RULISON  are  as  follows: 

(Gutenberg,  A>16‘’)  m^^  -  4.62  +  0.36  for  11  stations, 

(Gutenberg,  ^>15®)  -  3.99  +  0.44  for  12  stations. 
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IDENTIFICATION  CRITERIA  APPLIED  TO  RULISON 


Introduction 

One  of  the  main  objectives  of  the  VELA  program  has  been  to 
develop  seismological  criteria  for  distinguishing  between  under¬ 
ground  nuclear  explosions  and  earthquakes.  Prior  to  RULISON,  con¬ 
siderable  effort  had  been  devoted  to  establishing  diagnostic  crite¬ 
ria  making  use  primarily  of  Nevada  Test  Site  (NTS)  explosions, 
nearby  earthquakes,  and  earthquakes  from  both  Western  United  States 
and  elsewhere.  The  RULISON  explosion  provided  an  excellent  oppor¬ 
tunity  to  determine  the  general  applicability  of  these  diagnostic 
criteria  to  a  new  source  region.  A  chemical  explosion  was  detona¬ 
ted  in  the  same  general  region  as  RULISON  on  23  May  1967  at  Climax, 
Colorado;  however,  the  seismic  energy  released  was  not  sufficient 
to  be  detected  at  teleseismic  distances  thereby  seriously  limiting 
the  application  of  the  diagnostic  criteria  for  that  event.  Where- 
ever  possible,  these  criteria  have  been  applied  to  the  CLIMAX  explo¬ 
sion  and  to  a  selected  Colorado  earthquake  as  well  as  to  RULISON. 

The  earthquake  occurred  on  9  August  1967  and  was  located  by  USC5GS 
approximately  280  km  east  of  RULISON.  This  earthquake  was  selected 
on  the  following  basis  from  the  Preliminary  Determinations  of  Epi¬ 
centers  (USC8GS)  listing,  1)  magnitude  approximately  equal  to  that 
of  RULISON,  2)  a  shallow  event  and,  3)  location  as  near  to  RULISON 
epicenter  as  possible  and  fitting  the  above  requirements.  The  9 
August  1967  earthquake,  along  with  other  Denver  earthquakes  at  or 
near  the  Rocky  Mountain  Arsenal  (RMA)  well,  is  documented  and 
discussed  by  Major  and  Simon,  1968.  Further,  the  location  used 
for  this  study  was  determined  with  depth  restrained  to  5  km  and 
with  travel-times  from  43  WWSS  and  VELA  stations,  USC8GS  Earth¬ 
quake  Data  Report,  50-67,  (See  Appendix  5  for  the  event  spatial 
and  temporal  parameters.) 

Since  we  know  the  location,  origin  time,  and  size  of  RULISON 
it  is  important  to  show  the  degree  to  which  RULISON  revealed  itself 
seismologically .  Therefore,  in  this  section  we  discuss  the  follow¬ 
ing  diagnostic  criteria  as  applied  to  RULISON:  (1)  location, 

(2)  depth  of  focus,  (3)  complexity,  (4)  M^  versus  mj^,  (5)  energy 
relationships  among  phases,  (6)  radiation  patterns,  and  (7)  other 
criteria . 
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Location 


P-wave  arrival  times  were  read  from  all  available  VELA  obser¬ 
vatories  and  LRSM  stations  film  seismograms,  and  from  USC^GS  WWSSN 
film  chips.  Where  the  film  was  not  available,  USC§GS  arrival  times 
reported  in  the  "Earthquake  Data  Reports"  were  used  in  order  to 
obtain  a  network  of  stations  with  the  best  distance  range  and  azi¬ 
muth  aperture.  Location  according  to  Chiburis  (1966)  employing  the 
Herrin  1968  travel-time  tables  was  used  for  all  the  location  deter¬ 
minations  . 

RULISON,  when  located  using  66  stations  with  depth  free  to 
vary,  shifted  6  km  on  an  azimuth  of  31“  from  the  actual  location 
with  a  depth  estimate  of  41  km.  Using  a  subset  of  33  stations 
having  about  the  same  azimuthal  distribution,  the  event  was  mis- 
located  about  10  km  along  an  azimuth  of  202“  with  depth  estimated 
at  59  km. 

Results  of  the  cepstral  analysis  revealed  the  source  at  the 
expected  depth.  Accordingly,  RULISON  was  relocated  with  depth  re¬ 
strained  using  the  66  station  set,  which  produced  an  error  of  1  km 
along  an  azimuth  of  263“.  In  the  case  of  the  33  station  net,  the 
epicenter  shifted  24  km  along  an  azimuth  of  208“  from  the  true 
epicenter. 

The  earthquake  was  located  using  all  ten  available  arrival 
times  and  by  restraining  depth  to  the  USC§GS  determined  value  of 
5  km.  It  shifted  27  km  on  an  azimuth  10“  from  the  input  USC5GS 
location.  The  ten  stations  that  recorded  the  earthquake  also  re¬ 
corded  RULISON.  Using  this  ten  station  network,  RULISON  shifted 
16  km  on  an  azimuth  of  38“  from  its  true  location.  A  comparison 
of  location  results  is  presented  in  Table  III.  The  CLIMAX  explo¬ 
sion  was  not  located  because  there  were  no  teleseismic  data. 

Thus,  we  see  that  location  estimates  for  the  RULISON  event  are 
highly  variable  (1  km  to  24  km)  depending  on  the  recording  stations 
and  travel  times  used  and  the  source  depth  restrictions  placed  on 
the  computations. 

It  is  desirable  to  determine  the  degree  to  which  the 
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travel- times  derived  £rom  one  event  are  valid  for  use  in  locating 
another.  Accordingly,  travel  time  corrections  were  determined 
(Chiburis,  1968b)  and  the  events  were  relocated  with  the  network 
of  ten  common  stations.  Using  the  corrections  determined  from  the 
earthquake,  RULISON  shifted  15  km  on  an  azimuth  of  172°,  yielding 
no  improvement  in  accuracy.  We  would  expect  a  standard  deviation 
of  less  than  0.5  seconds  for  a  good  solution  with  corrected  travel 
times,  however,  the  standard  deviation  of  the  travel  time  residuals 
was  only  reduced  from  1.49  seconds  for  the  uncorrected  solution  to 
1.07  seconds  for  the  corrected  solution  with  a  corresponding  de“ 
crease  in  the  standard  95%  confidence  ellipse  from  4607  km^  to 
2396  km  .  Both  ellipses  comfortably  contain  the  actual  ejiicenter. 
The  fact  that  the  location  did  not  improve  is  not  surprising  when 
one  considers  that  the  two  events  are  separated  by  a  distance  of 
more  than  275  km  including  the  Rocky  Mountains.  Locating  the  earth¬ 
quake  with  travel  times  determined  from  RULISON  produced  a  shift 
to  14.42  km  on  an  azimuth  of  348.7°,  which  is  nearly  a  mirror 
image,  as  expected,  (Chiburis,  1969)  of  the  shift  obtained  for 
RULISON  with  the  earthquake  travel  time  corrections.  That  the 
standard  deviation  of  errors  did  not  decrease  significantly  demon¬ 
strates  that  that  the  travel  time  corrections  used  are  not  appro¬ 
priate.  The  conclusion  is  that  RULISON  and  the  Colorado  earthquake 
are  not  in  the  same  travel-time  region.  Although  the  area  in  the 
vicinity  of  RULISON  is  now  calibrated  for  locating  subsequent 
events,  the  area  in  the  vicinity  of  the  earthquake  remains  unknown 
regarding  suitable  corrections  for  accurate  locations. 

Depth  of  focus 

Since  the  depth  of  focus  for  explosions  is  presently  limited 
by  drilling  technology  to  shallow  depths  (less  than  10  km),  and 
since  most  earthquakes  are  known  to  occur  at  greater  depths  than 
this,  depth  determination  is  an  important  identification  criterion. 
Events  determined  to  be  significantly  deeper  than  a  few  kilometers 
can  be  classified  as  earthquakes.  Determination  of  the  focal  depth 
is  aided  by  the  fact  that  many  earthquakes  have  observable  depth 
phases  (pP,  sP,  etc.)j  explosions,  however,  are  generally  too 
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shallow  to  produce  visually  distinct  depth  phases  on  the  seismo¬ 
gram.  The  RULISON  event,  for  example  showed  no  depth  phase  on  the 
indivdual  recordings. 

For  shallow  events  the  pP  and  P  phases  overlap.  However,  if 
their  waveforms  are  similar,  the  fact  that  an  echo  exists  in  the  P 
coda  may  possibly  be  detected  by  the  presense  of  periodic  nulls  in 
the  P  spectra.  Cohen  (1969)  found  nulls  in  P  spectra  for  several 
explosions  due  to  the  interference  of  the  pP  and  P  phases.  By  com¬ 
puting  the  spectrum  of  the  spectrum  (the  cepstrum),  it  was  possible 
in  some  cases  to  estimate  objectively  the  depth-phase  delay  time. 
The  dot  product  of  the  cepstrum  and  the  pseudo-autocorrelation 
functions  shows  the  depth  phase  if  present  as  a  negative  correla¬ 
tion  peak. 

RULISON  and  the  9  August  1967  earthquake  are  two  shallow 
events  which  are  deep  enough  to  estimate  their  depths  by  the  spec¬ 
tral  null  method.  We  would  expect  RULISON  (2.6  km  depth)  to  gener¬ 
ate  a  pP  phase  1.5  seconds  20%  depending  upon  the  overburden 
velocity)  behind  the  P  phase,  with  spectral  nulls  every  0.7  Hz.  We 
would  expect  the  earthquake  at  S.O  km  depth  to  generate  a  pP  phase 
2.8  seconds  behind  the  P  phase,  and  with  spectral  nulls  every  0.4Hz. 

The  best  evidence  of  regular  spectral  nulls  and  of  the  promi¬ 
nent  negative  dot-product  correlation  peaks  are  in  the  LASA  data 
Figures  11a,  11b,  and  11c.  The  data  show  Pn  spectra  (0  to  4  Hz) 
for  RULISON  recorded  on  the  center  seismometers  of  18  LASA  sub¬ 
arrays.  The  dot  product  of  the  pseudo-autocorrelation  and  the 
cepstrum  accompanies  each  spectrum.  The  spectra  are  arranged  in 
order  of  increasing  distance  from  RULISON. 

The  sum  spectrum  (bottom.  Figure  illc)  gives  the  clearest  indi¬ 
cation  of  the  spectral  nulling  occurring  every  0.7  Hz  interval.  In 
addition,  the  negative  correlation  peak  at  1.45  seconds  is  the  most 
prominent  feature  on  the  dot  product  function.  These  results  agree 
with  the  expected  values  for  RULISON.  Many  of  the  individual 
spectra  show  some  of  the  spectral  nulls,  but  none  of  them  as  clearly 
as  the  summation  spectra.  In  addition,  many  of  the  individual  dot 
product  functions  show  the  negative  correlation  at  the  correct 
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delay  but  also  display  other  more  prominent  positive  and  negative 
correlation  peaks  not  associated  with  the  expected  depth  phases. 

LASA  data  for  the  9  August  1967  earthquake  is  not  available 
for  comparison  with  these  LASA  results  for  RULISON. 

Sum  spectra  from  LRSM  stations  for  both  RULISON  (17  stations 
at  Pn  and  P  distances)  and  the  9  August  1967  earthquake  (15  sta¬ 
tions  at  all  distances)  are  shown  in  Figure  18,  These  spectra  are 
log  plots  of  the  sums  shown  in  Figures  12c  and  13c.  Spectral  nulls, 
indicated  for  both  events,  are  more  prominent  for  the  earthquake 
spectrum.  Note,  too,  that  since  the  earthquake  is  deeper,  more 
spectral  nulls  are  observed  in  the  signal  band  of  0.5  Hz  to  4.0  Hz. 

Figures  12a,  12b,  and  12c  show  individual  and  summation  spec¬ 
tra  and  dot  product  traces  for  RULISON  obtained  from  recordings  at 
stations  at  all  distances.  The  regular  spectral  nulls  and  negative 
dot-product  correlation  peaks  are  decidedly  less  prominent  than 
those  shown  on  the  LASA  data. 

Figures  13a,  13b,  and  13c  show  individual  and  sum  spectra,  and 
dot  product  traces  from  the  9  August  1967  earthquake  for  15  LRSM 
stations  at  all  distances.  Again,  these  LRSM  stations  do  not  exhibit 
the  regular  spectral  nulls  and  prominent  negative  correlation  peaks 
as  well  as  does  the  LASA  data  for  RULISON. 

Conclusions 

1)  Regular  spectral  nulls  and  negative  correlation  peaks  were 
seen  for  both  events  where  expected  at  several  stations. 

2)  Not  all  stations  show  the  regular  spectral  nulls  and  cor¬ 
relation  peaks. 

3)  Nearly  all  stations  showed  other  spectral  nulls  and  nega¬ 
tive  correlation  peaks  besides  those  correlated  with  depth  phases. 

4)  LASA  data  exhibited  the  clearest  estimates  of  depth 
through  the  analysis  of  regular  spectral  nulls  and  dot-product 
correlation  peaks. 

5)  This  method  requires  summation  spectra  or  correlations 
from  several  traces  to  cancel  path  and  site  effects.  Multichannel 
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stations  like  LASA  are  helpful  in  this  way  since  the  structure 
varies  across  the  array. 

6)  RUL I S0\  showed  four  spectral  nulls  across  the  short-period 
signal  band  (0  5  Hz  to  3.0  Hz).  A  fewer  number  of  nulls  would 
make  recognition  of  regular  spectral  nulling  difficult.  Thus, 
events  with  depths  shallower  than  RULISON  would  be  difficult  to 
determine  by  any  variation  of  this  method. 

8)  The  9  August  1967  earthquake  showed  six  regular  nulls 
across  the  signal  band.  As  a  consequence,  the  depth  phase  for  this 
event  was  easier  to  detect  on  LRSM  station  data  than  was  RULISON* s. 

9)  The  method  requires  good  signal- to-noise  ratios. 
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Complexity 

The  generally  agreed  upon  requirements  for  the  satisfactory 
application  of  complexity  as  a  diagnostic,  (large  teleseismic  dis¬ 
tances,  good  signal  to  noise  ratios)  could  not  be  satisfied  because 
of  a  lack  of  stations.  However,  the  complexities  for  RULISON  and 
the  9  August  1967  earthquake  are  shown  in  Figure  14  plotted  as  a 
function  of  distance.  Both  events  have  similar  complexity  values 
(F^).  However,  for  common  stations;  LC,  KN,  SJ,  PG,  and  NP  the 
earthquake  has  greater  values  of  F^  except  at  NP  where  the  noise 
for  the  quake  is  so  high  that  one  can  hardly  see  the  signal. 
Basically,  the  signal  to  noise  ratio  is  inadequate  to  apply  this 
criterion. 

M  versus  m. 
s _ b 

Since  body  wave  and  surface  wave  magnitude  estimates  give  some 
measure  of  the  source  type  and/or  depths  by  minimizing  azimuthal 
and  distance  effects,  the  relative  excitation  of  body  and  surface- 
wave  energy  is  an  important  diagnostic  in  that  most  earthquakes 
have  larger  surface-wave  magnitudes  than  explosions  with  comparable 
body-wave  magnitude.  RULISON  by  virtue  of  its  new  location  and 
greater  depth  provided  an  opportunity  to  determine  the  effective¬ 
ness  of  this  criterion  as  compared  to  nuclear  explosions  from  the 
Nevada  Test  Site  (NTS)  and  earthquakes  from  the  Western  United 
States . 

Figure  15  shows  surface-wave  magnitudes  with  the  Gutenberg 
formula  applied  to  all  distances,  relative  to  Adjusted  body-wave 
magnitudes  (Adjusted  m^),  Evernden  (1967),  for  39  Nevada  Test  Site 
(NTS)  explosions.  For  comparison  we  also  include  twelve  earth¬ 
quakes  from  Western  United  States  with  magnitude  estimates  deter¬ 
mined  by  Basham.  (1969)  using  the  Canadian  station  network.  Since 
the  Canadian  seismological  network  includes  three  stations  at 
distances  less  than  15®  from  the  Nevada  Test  region,  we  excluded 
16  other  earthquakes  with  an  m^  determined  by  nine  or  less  stations 
for  the  purpose  of  minimizing  the  near  distant  station  effects. 
Therefore,  Basham's  m^^  should  be  comparable  to  the  Adjusted  m^^. 
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If  these  earthquakes  are  representative  of  th^  area  then  thlsy  fall 
in  a  different  population  than  the  explosions:  For  26  NTS  explo¬ 
sions  observed  teleseismically,  Basham  obtained  a  least  squares 
magnitude  relationship  of  =  1.24  m^  -  1.76.  We  obtain  M  =  1.21 
"'b  "  1*8  9  for  39  explosions  including  regional  and  near  regional 
observations.  Therefore,  the  latter  relationship  appears 

valid  for  NTS;  however,  the  separation  between  our  expiy)sion  and 
Basham's  earthquake  populations  is  enhanced.  This  is  largely  due 
to  the  differences  in  the  formulation  (0.18)  between  Basham  and 
Geotech  (see  Magnitude  section).  Also  included  in  the  figure,  but 
not  in  the  least  squarels  determination,  are  the  RULISON,  GASBUGGY 
and  CLIMAX  explosions,  the  23  January 'l966  New  Mexico  and  9  August 
1967  Colorado  earthquakes. 

As  shown  in  Figure  15,  M^  versus  Adjusted  mj^  the  three  explo- 
I  sions,  all  from  different  source  regions,  !fall  into  the  NTS  explo-  , 
sion  population  and  the  two  earthquakes  fall  in  the  group  of  other 
earthquakjes  in  Western  United  States.,  Thus  had  we  not  known  in 
advance  that  these  events  were;  explosions,  we  would  certainly  have 
decided  that  they  were  suspicions  on  the  basis  of  this  criterion 
alone. 

'  '  \  '  ' 

AR2  and  ERZ|  (area  ^nder  the  Rayleigh  wave  and  total  energy 

contained  in  the  Rayleigh  wave),  Turnbull  and  Lambert  (1968),  rel¬ 
ative  to  body-wave  magnitude  are  essentially  measiures  of  M^  1-  mj^. 
Figure  16  and'|17  show  ARZ  and  ERZ  plotted  versus  (Gutenberg)  mag¬ 
nitude  (mjj)  for  earthquakes  and  explosions  from  different  regions. 
The  tetm  (Gutenberg  mj^)  in  not  a  true  Gutenberg  magnitude  but 
includes  data  for  distances  less  than  16**  corrected  according  to 
VSC,  (see  Magnitude  section) .  These  figures  are  taken  from  the 
LONG  SHOT  report,  LaUert  et  al  (1970)  and  the  event  descriptions 
are  listed  in  Appendix  5.  The  results  place  RULISON  clearly  in 
the  (explosion  population.  1 

\ 

Thus,  on  the  basis  of  these  results  RULISON  would  fall  in  the 
general  explosion  population  as  defined  by  previous  NTS  (and  other) 
explosions  and' earthquakes^  This  does  not  necessarily  mean  that 
we  have  identified  RULISON  jas  an  explosion  at  this  point  but  rather 
that  jit  would  not  have  been  dismissed  as  being  an  earthquake. 
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Energy  relationships  among  phases  ^ 

\ 

The  previous  section  (M^  versus  shpws  that  a  difference 
in  excitation  of  body  and  surface  waves  exists  between  RULISON  and 
the  9  August  1967  Colorado  ea'rthquake.  The  purpose  of  this  section^ 
is  to  examine  these  differences  in  morp  detail.  Various  authors 
have  shown  that  more  long-period  energy  is  released  by  earthquakes 
than  by  explosions  having  equivalent  m^^;  however,  the  spectra  of 
earthquakes  seen  at  far-f ield ^are  influenced  by  source  orientation, 
depth,  and  distance  (see  Von  Seggern  and  Lambert,  1969,  for  the 
bibliography  and  summary  discussion).  Therefore,  we,  discuss  (1)  \ 

P-wave  spectra,  (2)  I^ayleigh-wave  spectra,  (3)  spectral  energy 
ratios,  and  (4)  relative  shear  energy. 

(1)  P-wave  spectra  - 

\ 

Figure  18,  discussed  in  the  Depth  of  Focus  section  ijs  a  log 
plot  of  the  sum  of  17  stations  for  RULISON  and  15  stations  for 
the  9  August  1967  earthquake  shown  in  Figuris  12c  i^nd  13c  respec¬ 
tively.  From  this  display  the  earthquake  does  show  more  long- 
period  energy  than  the  explosion.  This  is  in  qualitative  agree- 
with  most  resvilts  of  other  works  in  the  field  of  short-period  i 
spectral  ratios. 

I 

,  i 

(2)  Rayleigh  wave-spectra  - 

Rayleigh-wave  energy  spectra  are  computed  by  a  Isine-cpsine 
Fourier  transform'  for  a  signal  velocity  window  of  about  3.60  km/sec 
to  2.5  km/sec  and  a  noise  sample  of  equal  length.  The  noise  spec¬ 
trum  ip  subtracted  from  the  signal  spectrum.  The  resultant  i^nergy 
spectrum  is  then  reduced  to  an  amplitud^  spectrum  which  is  further 
corrected  for  static  magnification  and  system  response.  Figures 
19a  and  19b  show  Rayleigh-wave  spectra  of  five  commlon  stations 
(KN-UT,  LC-NM,  SJ-TX,  RK-ON,  and  NP-NT)  for  RULISON  and  9  August 
1^67  earthquake.  Fourteen  additional  RULISON  spectra  are  shown  in 
Appendix  6. 

In  Figure  20  We  show  summed  and  normalized  energy  spectra  for 
these  same  five  common  stations  for  each  event,  ihis  process  of 
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summing  Rayleigh-wave  spectra  smears  out  the  spectrum  because  of 
frequency  dependent  attenuation  with  distance.  Since  the  RULISON 
and  earthquake  epicenters  are  separated  by  a  large  distance 
(280  km)  this  method  could  also  introduce  large  path  effects.  How 
ever,  it  is  interesting  to  note  that  the  spectral  sums  are  very 
similar  in  shape  between  0.02  and  0.076  cps.  The  earthquake  does 
show  more  energy  in  the  low  (0.02  to  0.38  cps)  and  high  (0.076  to 
.10  cps)  frequency  portions,  however.  This  result  is  discussed 
more  quantatively  in  the  next  section. 


(3)  Energy  ratios  - 

Von  Seggern  and  Lambert  (1969)  studied  the  Rayleigh  wave 
spectral  dependency  as  functions  of  magnitude,  source  type,  and 
distance  both  theoretically  and  empirically.  The  theory  indicates, 
for  explosions,  that  the  shape  of  the  source  spectra  should  not 
change  with  magnitude  in  the  spectral  band  of  interest  (T  -  10  to 
50  seconds);  thus,  ratios  of  total  energy  between  adjacent  bands 
of  frequency  should  not  change  with  magnitude.  However,  due  to 
different  source  types,  layering  response  and  depths,  earthquake 
ratios  can  be  greater  than,  less  than,  or  equal  to  those  of  explo¬ 
sions.  Empirical  evidence  supports  these  theoretical  findings. 

Mean  ratios  for  RULISON,  the  9  August  1967  Colorado  earth¬ 
quake  and  the  23  January  1966  New  Mexico  earthquake  are  computed 

as  follows: 

’^1  -  K  ?  I  t 

m=l  '  t  1 2 

where  1-  equals  the  average  energy  ratio  for  n  stations  recording 
one  e^^ent  and  is  the  energy  spectra  at  station  m.  For  we  let 

f^  =  1/T^  and  T^  =  48  seconds, 

^2  “  ^^^2  ^2  "  seconds, 

fj  =  l/Tj  and  Tj  =  15  seconds. 

For  ^2,  we  change  only  T3  from  15  to  10  seconds  (T3  =  10  seconds). 
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Figures  21  and  22  show  and  ^2  plotted  versus  Rayleigh-wave  mag¬ 
nitude  and  the  event  descriptions  are  listed  in  Appendix  7.  The 
Rayleigh-wave  magnitude  is  determined  from  the  spectra  used  in  the 
above  energy  ratio  analysis.  For  these  estimates  the  maximum  am¬ 
plitude  is  picked  from  the  spectra  between  T  =  17  and  23  seconds 
and  corrected  according  to  Gutenberg  formulation  applied  to  all 
distances.  Therefore  these  magnitudes  do  not  correspond  to  those 
determined  from  film  analysis.  It  should  be  emphasized  that, 
basically,  the  ratios  are  independent  of  magnitude,  but  only  if 
the  signal  to  noise  ratio  is  high  enough  to  allow  a  valid  analysis. 

In  Figure  21,  R^  (T^  =  15  sec)  determined  from  19  stations  for 
RULISON  is  clearly  in  the  explosion  population.  Similar  results 
for  R2  (T^  =  10  sec)  are  obtained  in  Figure  22. 

determined  from  14  stations  for  the  Colorado  earthquake  is 
greater  than  for  RULISON  and  places  it  in  the  earthquake-popula¬ 
tion  while  ^2  foi^  fhe  earthquake  places  it  in  the  mixed  explosion- 
earthquake  population.  Similar  results  using  14  stations  were 
obtained  for  the  New  Mexico  earthquake.  However,  Von  Seggern  and 
Lambert  (1969)  indicate  that  for  events  with  greater  than  3.0 
and  an  average  station  network  distance  greater  than  1000  km,  R^ 

(T3  “  15  sec)  is.  more  reliable  as  a  diagnostic  than  T?2  (T^  =  10 
sec).  The  average  station  network  distances  for  these  three  events 
ranged  from  1700  to  2000  km. 

Therefore,  for  RULISON,  the  preferred  diagnostic,  shows 
that  the  Rayleigh  wave  spectra  in  the  T  =  15  to  50  second  band  is 
similar  to  other  explosions. 

(4)  Relative  shear  energy  - 

Earthquake  source  mechanisms  by  virture  of  their  physical  con¬ 
figuration  should  produce  more  shear  energy  than  compressive  source 
mechanisms  such  as  explosions.  Thus,  we  look  at  relative  Lg  and 
Love  wave  amplitudes  for  the  9  August  1967  earthquake  and  RULISON. 

The  average  ratio  of  Lg  earthquake  amplitudes  to  Lg  RULISON 
amplitudes  is  6.62  and  the  average  ratio  of  Love  earthquake 
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amplitudes  to  Love  RULISON  amplitudes  is  6.58. 

Therefore  the  earthquake  has  more  shear  energy  than  RULISON, 
an  equivalent  magnitude  explosion.  However,  additional  data  from 
many  Colorado  earthquakes  will  be  needed  to  ascertain  whether  this 
is  typical  or  atypical  of  this  region.  Moreover,  there  is  no  well 
accepted  discriminant  function  for  this  ratio  against  which  to  test 
these  two  data  points. 

Radiation  patterns 

Various  authors  have  shown  that  differences  in  radiation  pat¬ 
terns  between  earthquakes  and  explosions  should  exist.  Since  path 
and  site  effects  greatly  influence  body-wave  amplitudes  such  that 
definitive  radiation  patterns  are  rarely  possible,  we  look  at  only 
the  direction  of  first  motion  for  P. 

For  RULISON  and  the  Colorado  earthquake  (9  August  1967)  there 
appeared  to  be  no  distinct  rarefactional  first  motion;  however,  the 
quality  of  the  seismograms  prevents  any  reliable  determination  of 
first  motion.  Thus,  neither  the  earthquake  nor  RULISON,  could  be 
dismissed  as  an  earthquake. 

Studies  by  Toksoz  et  al  (1965)  and  Toksoz  and  Clermont  (1967) 
have  shown  that  observations  of  non-circular  radiation  patterns  of 
surface  waves  for  the  nuclear  explosions,  HARDHAT,  HAYMAKER,  SHOAL, 
and  BILBY  can  be  explained  in  terms  of  a  compressional  source  accom¬ 
panied  by  the  release  of  tectonic  strain. 

For  the  theoretical  aspects  and  procedural  considerations  of 
this  study  we  reference  Toksoz  and  Clermont  (1967)  and  show  formula 
7  in  the  above  paper. 


U, 


Fk^Aj^  cos2e 


|^Rz|  (1+F  sin  20)Aj^k^(uJ/w^) 


exp  [-r(YL  -  Yr)  ] 


displacement  of  Love  waves, 

far-field  displacement  of  Rayleigh  waves, 

relative  strength  of  the  double-couple, 

azimuthal  dependence  of  the  Love  wave  radiation  for 
the  double-couple, 

azimuthal  dependence  of  the  Rayleigh  wave  radiation 
for  the  double-couple, 

wave  numbers  for  Love  and  Rayleigh  waves, 
radial  distance. 

Love  and  Rayleigh  wave  attenuation  coefficients, 

the  medium  response  for  Love  and  Rayleigh  waves  due 
to  a  vertical  force, 

the  components  of  particle  velocity  at  the  surface. 

We  determined  Aj^,  Aj^,  kj^,  kj^,  and  Cu*/w*)  using  Harkrider's 
program  (1964)  for  an  average  structure  of  the  Western  United 
States  as  given  by  Alexander  (1963).  These  parameters  were  deter¬ 
mined  for  a  period  (T)  of  13.0  seconds  since  most  of  the  observed 
peak  amplitudes  of  the  Love  and  Rayleigh  waves  for  RULISON  were 
near  this  value  (Table  II). 

Figure  23  shows  all  possible  data  of  U^^/Uj^  for  RULISON  and 
also  the  theoretical  radiation  pattern  based  on  a  strike  plane  of 
the  double-couple  system,  0  =  24°,  and  a  relative  strength  (F)  of 
the  double-couple  to  explosion  of  0.60. 

The  strike  angle  is  normal  to  the  strike  of  the  tectonic  fea¬ 
tures  in  the  vicinity  of  the  RULISON  epicenter  (Figure  24).  Toksoz 
indicates  that  the  relative  strength  (F)  seemingly  depends  upon  the 
properties  of  the  source  medium  and  cites  the  following: 

(1)  GNOME  and  SALMON  (salt),  SEDAN  (loose  aluvium);  f  =  0.0. 

(2)  HAYMAKER  (allivium);  F  =  0.333. 

(3)  BILBY  (tuff);  F  =  0.47. 

(4)  SHOAL  at  a  relatively  greater  depth  in  granite;  f  =  0.90. 
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On  this  basis  alone  the  F  “  0.60  obtained  for  RULISON,  perhaps 
appears  greater  than  expected.  It  is  possible  that  depth  and  the 
immediate  stress  field  as  well  as  medium  would  influence  this  value. 

A  further  consideration  of  importance  to  the  identification 
problem  is  a  comparison  of  Love  to  Rayleigh  wave  amplitude  ratios 
between  RULISON  and  9  August  1967  Colorado  earthquake.  Figure  25 
shows  the  LQ/LR  amplitude  ratios  for  the  Colorado  earthquake.  Even 
though  we  did  not  determine  the  earthquake  double-couple  parameters, 
it  is  evident  that  the  differences  are  dramatic,  as  compared  to 
RULISON  both  in  radiation  pattern  and  magnitude  of  ratios.  The 
average  LQ/LR  amplitude  ratio  for  RULISON  is  0.52  and  for  the 
earthquake  3.48.  Major  and  Simon  (1968)  indicate  right  lateral 
fault  movement  determined  from  RMA  close-in  strainmeter  records  for 
the  9  August  1967  earthquake.  Further,  the  event  originated  at  a 
depth  of  5.0  km  about  2.7  km  north  and  4.6  km  west  of  the  RMA  well 
and  propagated  N70®W  for  a  distance  of  about  10  km  with  a  velocity 
of  about  3.0  km/sec. 

Thus,  it  is  clear  that  the  two  events  have  distinctly  differ¬ 
ent  source  mechanisms;  however,  additional  data  from  many  Colorado 
earthquakes  will  be  needed  to  ascertain  whether  this  earthquake  is 
typical  or  atypical  of  this  region. 

•II 

Other  criteria 

Seismic  events  detected  from  aseismic  regions  and  far  from 
known  nuclear  test  areas,  are  ones  that  would  naturally  be  subjec¬ 
ted  to  further  study  for  the  purpose  of  identification.  For  such 
events  the  presence  or  absence  of  aftershock  activity  would  also 
be  of  considerable  importance. 

RULISON  occurred  in  northern  Colorado  away  from  a  known  test 
region.  The  USC8GS,  PDE  listing  was  searched  for  events  occuring 
in  a  2°  (latitude  and  longitude)  area  around  the  RULISON  epicenter 
between  9  August  1966  and  1  January  1970.  No  seismic  events  were 
reported  except  RULISON  (Figure  26).  Further,  from  this  listing 
and  film  analysis,  no  aftershocks  were  observed.  In  addition,  the 
origin  time,  as  determined  from  the  location  results.  Table  III, 
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is  within  2.5  seconds  of  the  21:00:00.02  hour.  In  contrast,  a 
similar  2®  area  was  studied  for  the  9  August  1967  earthquake  and 
14  events  were  reported.  However,  all  but  one  of  these  events 
should  also  be  considered  suspicious  because  of  the  reported  shal¬ 
low  depths  (5  km)  and  close  spatial  grouping,  near  the  RMA  well. 
Table  IV  -and  Figure  27. 

Thus,  RULISON  is  located  in  an-  aseismic  uncommon  test  region 
with. an  origin  time  of  21:00:00.02. 

Summary  and  conclusions 

The  following  table  summarizes  the  results  of  applying  the 
various  discriminants  discussed  above  to  RULISON.  i 


SUMMARY  OF  RULISON  DIAGNOSTICS 
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Figure  1.  Station  location  map  for  RULISON. 
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Figure  2.  Reduced  travel -times  for  RULISON 
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Figure  6.  LQ  amplitudes  for  RULISON. 
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Figure  7.  LR  amplitudes  for  RULISON. 
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Figure  8,  Body-wave  magnitudes  for  RULISON. 
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Figure  10.  Short  pe 
prof i 1e. 
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Figure  Ha.  Cepstrum  analysis  for  RULISON  at  LASA  sub-array 
centers.  ^ 
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Figure  lib.  Cepstrum  analysis  for  RULISON  at  LASA  sub-arrav 
centers. 
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Figure  11c.  Cepstrum  analysis  for  RULISON  at  LASA  sub-array 
centers . 
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VERTICAL  SCALE  VALUES  NORMALIZED  TO  PLOT  SCALE  OF  1  INCH 
Figure  12b.  Cepstrum  analysis  of  Pn  and  P  for  RULISON. 
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Figure  13a.  Cepstrum  analysis  of  Pn  and  P  for  the  Colorado 
earthquake.'  - 
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Figure  13b.  Cepstrum  analysis  of  Pn  and  P  for  the  Colorado 
earthquake. 
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Figure  13c,  Cepstrum  analysis  of  Pn  and  P  for  the  Colorado 
earthquake. 
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Figure  14.  RULISON  complexity  versus  distance 
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Figure  15.  Me  (Gutenberg)  versus  1115  for  39  NTS  explosions  and 
14  western  United  States  explosions. 
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Figure  18.  Comparative  SP  spectral  energy  summation  for 
RULISON  and  9  August  1967  earthquake. 
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Figure  19*.  Coaiperetive  LR  inplltude  spectre  for  RULISON 
and  9  August  1967  eerthguaice  at  KN-UT,  LC-NN, 
SJ-TX.  RK-ON.  and  NP-NT. 
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Figure  19b,  Conperetive  LR  emplltude  spectre  for  RULISON 
end  9  August  1967  eerthqueke  et  KN>UT,  LC-NM, 
SJ-TX,  RK-ON.  end  NP-NT. 
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Figure  20.  spectre)  energy  lUMiatlon  for  RULISON 

end  9  August  1967  eerthqueke. 
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Figure  24,  Map  of  the  Colorado  basement  features. 
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Figure  26.  Seismicity  of  the  RULISON  source  region. 
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TABLE  I 

# 

Event  Description 
For 
RULISON 

Date:  10  September  1969 

Time  of  Origin:  21:00:00.12 

Magnitude:  4.62  with  standard  deviution  of  +0.36  (11  stations) 

Location:  Western  Colorado  situated  between  the  White  River 

uplift  and  the  Uinta  Basin 

Coordinates : 

Latitude:  39‘’24'2r'N 

Longitude:  107®56'53"W 

Environment: 

Geological  Medium:  Mesa  Verde  Sandstone 
Surface  Elevation:  8154  feet  above  mean  sea  level 
Shot  Elevation:  277  feet  below  mean  sea  level 

Shot  Depth:  8431  feet  below  surface 
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TASU  11  (Cant'd.) 
Arlttclpnl  Pfetaaa  far  DULISON 


coot  STATION  (t;r  iS.,,. 
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List  of  Events  Occurring  in  the  9  August  1967  Colorado  Earthquake  Source  Region 
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APPENDIX  I 

Station  Locations  and  Elevations  for  RULISON 
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